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Abstract. Some direct detection experiments have recently collected excess events that could 
be interpreted as a dark matter (DM) signal, pointing to particles in the ~10 GeV mass range. 
We show that scenarios in which DM can self-annihilate with significant couplings to quarks 
are likely excluded by the cosmic-ray (CR) antiproton data, provided the annihilation is S- 
wave dominated when DM decouples in the early universe. These limits apply to most of 
supersymmetric candidates, e.g. in the minimal supersymmetric standard model (MSSM) and in 
the next-to-MSSM (NMSSM), and more generally to any thermal DM particle with hadronizing 
annihilation final states. 



1. Introduction 

After the finding of an annual modulation consistent with a DM signal in the DAMA/Libra 
data [1, 2], some skepticism had arisen due to the failure of other direct detection experiments 
to confirm it in the preferred phase space corner. Interestingly enough, this initial caution has 
been recently tickled by the results of a few other experiments with different setups. CoGeNT ^ 
and CRESST ^ have indeed found excess events; CoGeNT has even released hints for an annual 
modulation [5]. Yet, XENON-10/100 [61 [7], CDMS [8] and SIMPLE [9] do not validate such 
excesses and therefore challenge the classical DM interpretation of the aforementioned positive 
results — i.e. the elastic collisions of weakly massive interacting particles (WIMPs) off target 
nuclei [10]. It is still important to recall that these excesses — and bounds — rely on data 
analyses led very close to the current experimental thresholds, where systematic errors can in 
principle get large and remain to be much better controlled. Nevertheless, these measurements 
have brought some excitement in the DM community, and many authors have tried to fit the 
available data with some phenomenological models {e.g. [H]). The preferred mass range is 
around 10 GeV, while difficulties are generally found when trying to get good fits to the signals 
while avoiding the constraints from other direct detection experiments |10j . Some ideas have 
been proposed to explain such discrepancies {e.g. [I2l [l3]), but what people usually do is to 
forget about either the bounds or the signals, depending on whether they want to promote their 
models or challenge them. 
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Figure 1. Left: Impact of considering 2 different temperatures for the quark-hadron transition 
on the relic density calculation as a function of the WIMP mass, for a given annihilation 
cross section [16]. Middle: antiproton flux prediction for a 10 GeV WIMP annihilating into 
bb [15]. Right: antiproton flux predictions for a 12 GeV WIMP annihilating into different mass 
combinations of an intermediate two-boson state which further decays into quarks [16j. 



Beside direct detection, indirect detection, that looks for DM annihilation or decay traces 
originating in high concentration astrophysical sites, provides an interesting way to check 
the relevance of any phenomenological proposal in this context. Among indirect signatures, 
CR antiprotons have long been known as an interesting target to look at jT3j. In this 
proceedings' contribution, we show that they indeed provide quite strong limits on quarkophilic 
DM candidates, and are even likely excluding those models which have an S-wave dominated 
annihilation cross section. The following discussion is mostly based upon two articles, Refs. |15j 
and ^6], to which we refer the reader for more details and references. 

2. Ingredients and results 

Before discussing the Galactic ingredients, it might be useful to recall a source of uncertainty 
in the annihilation cross section that comes from the calculation of the relic density, though 
commonly assumed to be under control. For typical thermal cold DM candidates, the 
chemical decoupling that sets the DM abundance occurs in the early universe at a temperature 
Tf ~ m^/20. In the 10 GeV mass range, this translates into temperatures around 500 MeV, 
a temperature range in which the quark-hadron transition is expected [17| 118] . A decoupling 
occurring before or after this transition induces a factor of 2 uncertainty in the annihilation 
cross section for a given cosmological abundance, since the relativistic degrees of freedom, which 
characterize the energy density in the radiation era (i.e. the squared expansion rate), vary by a 
factor of 4 meanwhile. We therefore note as a preliminary remark that the required annihilation 
cross section may be twice larger than the canonical value in this mass range, the transition 
mass being as uncertain as the quark-hadron transition temperature. This was discussed in 
detail in |16j , and is illustrated in the left panel of Fig. [T] We now return to the prediction of the 
CR antiproton flux, for which the involved physics is a little more complex than for gamma-ray 
calculations. The additional steps come from that charged CRs scatter off magnetic turbulences 
in the Galaxy (which drives their motion to a diffusion pattern), are convected away from the 
disk by Galactic winds, and suffer nuclear reactions (spallation) with the interstellar medium 
(ISM) as well as energy redistribution due to scattering interactions and diffusive reacceleration. 
For a complete theoretical view on CR propagation, see Ref. [19]. Thus, a full prediction relies 
on solving a diffusion equation that reads (steady state and stable species): 



-V{K^VM -VcJ^} + dE{jCEM^ + rsM = Q{x,E) 



(1) 



In this equation, N = dn/dE is the CR density per unit energy, Kx is the spatial diffusion 
coefficient, Vc is the convective wind velocity, is the rate of spallation with the ISM gas and 
KLe is an energy redistribution (differential) operator. The rhs term is the source term. For DM 
annihilation, it is proportional to the product of the annihilation cross section and the squared 
DM number density (the squared DM mass density to the squared WIMP mass), {av) (p/m^)^. 
Semi- analytic solutions to the previous equation do exist in certain regimes, especially when 
one fixes the geometry of the diffusion zone to a slab inside which the diffusion coefficient is set 
homogeneous. Further reasonably fixing the spatial extent of this slab, i.e. its radius i? ~ 20 
kpc and half-width L ~ 4 kpc, closes the full definition of the propagation model. All these 
parameters are globally referred to as transport parameters. They are usually constrained with 
CR nuclei data — basically, one predicts the fiuxes of those secondary species arising from 
the spallation of some primary astrophycical CR nuclei, the source of which may be rather 
well modeled (the a priori unknown source parameters cancel out at zeroth order when using 
the secondary-to-primary ratios — e.g. |20j). Predictions of the CR antiproton flux due to 
DM annihilation have been derived for years {e.g. [H]), but were often inconclusive because 
of uncertainties in the diffusion halo size, a diffusion zone down to L ~ 1 kpc inducing an 
unobservable flux. Nevertheless, the constraints on the diffusion volume have been strongly 
improved since then, and it is now expected to extend to L > 3 kpc [22j — it is also interesting 
to note that taking L I kpc inevitably leads to a strong excess of secondary astrophysical 
positrons below 5 GeV [23l 123] . The results to be presented below were essentially derived using 
L = 4 kpc, and a diffusion coefficient normalization of Kq = Kx{p/q = IGV) ~ 0.01 kpc^/Myr. 

Predictions of the DM-induced antiproton fiux are poorly sensitive to the distribution of DM 
close to the Galactic center — as long as p oc r~"' with 7^1-5 — in contrast to gamma-ray 
predictions. This is due to diffusion itself, that tends to smear out far away inhomogeneities. 
Actually, the antiproton fiux is much better controlled by the local DM density pq. There has 
been some recent attempts to bracket it, and it is noteworthy that quite different methods have 
led to similar results, pq ~ 0.4, with a ~ 20% accuracy {e.g. [25l [26l [27] ) . The coming results 
can be derived by using the best-fit Einasto or Navarro-Frenk- White PU] profiles provided 
in Refs. [25J and [27j, indifferently. 

As to DM models, we are going to consider two template cases that are likely representative 
of almost all quarkophilic configurations. In the first one [TS], annihilation fully goes into 
56-quark pairs, which is conservatively representative of annihilation into any quark pair, up 
to an antiproton multiplicity factor typically < 2 in the interesting ~ 0.1 — 5 GeV kinetic 
energy range, unfavorable to heavy quarks (smaller dynamical range). In the second one [16j, 
annihilation goes into two unstable particles, scalar and pseudo-scalar, with any masses above 
1 GeV allowed by kinematics, and which may further decay into quarks. These (pseudo-)scalar 
particles could be thought of as the light Higgs bosons arising in the NMSSM, or more generally 
in any singlet extension of the MSSM. Except for spin considerations, this annihilation channel 
is kinematically representative of any two-body intermediate state made of unstable bosons. The 
following analysis can therefore be extrapolated, for instance, to models that involve couplings 
to new light gauge fields [30] . 

The antiproton constraints merely come from the fact that the expected secondary 
astrophysical background already saturates the existing data, e.g. from BESS [31\ or 
PAMELA [32j. Consequently, any other significant astrophysical yield quickly leads to an 
excess prediction. The important point is that the theoretical errors affecting the astrophysical 
background estimates are small, at the 15% level below a few GeV [33]- It turns out that 
the DM-induced antiproton fiux generically exceeds the background expectation for WIMP 
masses between ~5-15 GeV and (av) > 10^^^ cm^/s, for the two configurations discussed above. 
Therefore, if the annihilation cross setion is S-wave dominated at the decoupling time in the 
early universe, then requiring the correct DM relic abundance, which translates into having 



2.5 X 10^^^ < {av) 5 X 10^^^ cm^/s, generically leads to an antiproton excess in this mass 
range — a smaller cross section is difficult to achieve, since leading to an over-abundance of 
DM. Examples are given in the middle and right panels of Fig. [T] The only ways left to escape 
these antiproton bounds are (i) to consider the second configuration, i.e. the annihilation into a 
two-boson intermediate state where the decaying bosons would both have masses lower than the 
(anti)proton mass, or (ii) to be at the 6-quark pair threshold production at which antiprotons 
can also barely be produced. The first case might be easy to cook up, but the second one is 
clearly very fine-tuned. 



3. Conclusion 

We have discussed how CR antiprotons can put limits on light DM candidates. They are 
generically strong and exclude quarkophilic DM particles in the 5-15 GeV mass range relevant 
to the direct detection signals, if they annihilate through S-waves with {av) > 10-26 cmVs. 
The only constraint at the origin of this result is the upper bound on the cosmological DM 
abundance, which translates into a lower limit on the annihilation cross section at the decoupling 
time, (av) > 2.5 x 10"^^ cm^/s. Light neutralinos arising in the MSSM with non-unified gaugino 
masses and in the NMSSM are therefore excluded as potential DM candidates if they fall in the 
quarkophilic categories discussed above, which is quite often the case. These antiproton limits 
are complementary to other studies using e.g. gamma-rays |3l], CMB photons [351, and solar 
neutrinos [36], which concern more annihilation final states than only quarks. 
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